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Abstract 

Several wildlife sanctuaries in the world are home to the surviving populations of many endemic 
species. Trishna wildlife sanctuary in northeast India is protected by law, and is home to the last 
surviving populations of Asian bison {Bos gorus Smith), spectacle monkey {Trachypithecus 
phayrie Blyth), capped langur {Trachypithecus pileatus Blyth), slow loris {Nycticebus coucang 
Boddaert), wild cat {Felis chaus Schreber), and wild boars {Sus scrofa L.), among many other 
animals and plants. The sanctuary was explored for species richness and diversity of butterflies. 
A six-month-long study revealed the occurrence of 59 butterfly species that included 21 unique 
species and 9 species listed in the threatened category. The mixed moist deciduous mature forest 
of the sanctuary harbored greater species richness and species diversity (39 species under 31 gen- 
era) than other parts of the sanctuary, which is comprised of regenerated secondary mixed 
deciduous forest (37 species under 32 genera), degraded forests (32 species under 28 genera), and 
open grassland with patches of plantations and artificial lakes (24 species under 17 genera). The 
majority of these species showed a distribution range throughout the Indo-Malayan region and 
Australasia tropics, and eight species were distributed in the eastern parts of South Asia, includ- 
ing one species, Labadea martha (F.), which is distributed in the eastern Himalayas alone. 
Estimator Chao 2 provided the best-predicted value of species richness. The steep slope of the 
species accumulation curve suggested the occurrence of a large number of rare species, and a 
prolonged gentle slope suggested a higher species richness at a higher sample abundance. The 
species composition of vegetation-rich habitats showed high similarity in comparison to vegeta- 
tion-poor habitats. 
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Introduction 



Recent studies of biodiversity in relation to 
ecosystem functioning have suggested that 
species diversity sometimes enhances produc- 
tivity and stabihty of ecosystems (Naeem et 
al. 1994; Tilman et al. 1996). Positive rela- 
tionships have been found between butterfly 
diversity and plant diversity (Thomas and 
Malorie 1985; Leps and Spitzer 1990). This 
relationship is particularly true in tropical re- 
gions, where insects show high abundance and 
species diversity (Price 1997). Due to richness 
in vegetation, the northeastern region of India 
is home to a rich diversity of butterflies, 
among other insects (Kunte 1997; Alfred et al. 
2002), and it is also part of one of the mega 
biodiversity hotspots of the world (Myers et 
al. 2000). A review of literature suggests that 
76 species of butterfly were previously rec- 
orded from the state of Tripura (10,490 sq 
km), northeast India (Mandal et al. 2002; 
Agarwala et al. 2010; Majumder et al. 2011). 
Among other northeastern states, 104 species 
of butterfly from Meghalaya (22429 sq. km), 
695 species from Sikkim (7096 sq. km), and 
962 species from Assam (78438 sq. km) have 
been recorded (Evans 1932; Talbot 1939; 
Wynter-Blyth 1957; Haribal 1992). Evidently, 
the knowledge base of lepidopteran fauna and 
their distribution in different habitats is une- 
ven and still scanty from this part of India. 

In this study, a detailed inventory was carried 
out to document the butterfly species richness 
and diversity of Trishna wildlife sanctuary 
(TWLS) in Tripura in relation to its habitat, 
which is comprised of mixed moist deciduous 
and evergreen forests with patches of long 
grasses. At the time of this study, 22 species 
were known from this sanctuary (Roy 
Choudhury et al. 201 1) without regard to their 
abundance, distribution pattern, and habitat 
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preference. Considering the fact that Trishna 
sanctuary is known to have 230 tree species, 
110 species of shrubs, 400 species of herbs, 
and 150 species of climbers ( www.incredible- 
northeastindia.com ), the assumption was that 
new data on spatial scale over a period of time 
would help to ascertain the true status of spe- 
cies richness and diversity in diverse habitats 
of the Trishna wildlife sanctuary, which is 
home to the Asian population of bison {Bos 
gorus Smith), and 116 species of resident and 
migratory birds, reptiles like cobra and king 
cobra snakes, and mammals like wild boar 
{Sus scrofa L.), spectacle monkey 
(Trachypithecus phayrie Blyth), capped lan- 
gur {Trachypithecus pileatus Blyth), slow 
loris {Nycticebus coucang Boddaert) and wild 
cat (Felis chaus Schreber) 

( www, tripura.nic . in/trishna/intro .html ) . This 
study also intended to bring out any hitherto 
not recorded threatened taxa of butterfly from 
this natural preserve of wildlife in this part of 
south Asia. 

Materials and Methods 

Study area 

Trishna wild life sanctuary is situated in the 
south district of Tripura state between 23"^ 
26.137' N and 91^ 28.184' E and has an alti- 
tudinal gradient of 51-82 m a.s.l. (Figure 1). 
The total sanctuary area is 194.71 km^and is 
delimited on the east and west sides by the 
international boundary with Bangladesh. The 
forest cover of the sanctuary consists of dense 
primary forest (62%) dominated by Shorea 
robusta Roth, Dipterocarpus turbinatus 
C.F.Gaertn, and Terminalia belliraca (Gaertn) 
trees, degraded forest (18%) dominated by 
Toona ciliata M. J. Roem, Albizia procera 
Durazz, a large number of shrubs, herbs and 
climbers, and the remaining 20% is bushy for- 
est of bamboos, sedges, long grasses, and 
shrubs like Microcos paniculata L. of Chinese 



Journal of Insect Science | http://www.insectscience.org 



2 



Journal of Insect Science: Vol. 1 3 | Article 79 

origin, Chromolaena odorata (L.) King and 
H.E. Robins of North American origin, and 
Lantana camara L. of Neotropical origin, 
among others. The sanctuary has a number of 
perennial water rivulets and patches of grass- 
lands. The climatic condition is tropical, with 
a minimum rainfall of 3.58 mm in December 
and a maximum of 508.20 mm in July. The 
average minimum and maximum temperature 
recorded in the region is 6.8"^ C in January and 
37.70"^ C in June, respectively. 

This study was conducted in four distinct hab- 
itats of TWLS based on habitat type, 
topography, and intensity of anthropogenic 
pressure. A two-day scouting survey was car- 
ried out in April 2010 to identify the sampling 
sites, one each in four habitat types (Table 1). 
TWLS I (near an ecotourism site) was charac- 
terized by mature, secondary mixed moist 
deciduous forest dominated by Shorea, Dip- 
terocarpus, and Microcos plants, with 
multilayered understory vegetation (sedges, 
grasses, and ferns) and a canopy height of 15 
m to 20 m. The major human pressures were 
fodder and fuel wood collection by inhabitants 
settled on the periphery of the sanctuary. Two 
artificial lakes were present on each side of 
TWLS L TWLS II (near a medicinal plant 
nursery) was characterized by regenerated, 
secondary mixed moist deciduous forest dom- 
inated by Dipterocarpus but few Shorea trees, 
and little understory, delimited by profuse 
growth of herbs and shrubs at the edges. The 
canopy height was between 3 m and 8 m, and 
the site was surrounded by paddy {Oryza sati- 
va L.) fields. Wild boars inhabited this site 
and caused uprooting of understory vegeta- 
tion. TWLS III (near a bison watch tower) 
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was characterized by degraded forest with 
growths of herbs, creeper grasses, and small to 
large patches of several bamboo species 
(B ambus a tulda Roxbergii, B. balcooa: Rox- 
bergii, and Meloccana baccifera: (Roxburgh) 
Kurtz ex Skeels) with very little understory 
vegetation. The forest floor was covered with 
leaf litter of 10 cm to 25 cm thickness. TWLS 
IV (a bison feeding ground) was devoid of big 
and old trees, and was an open type forest 
dominated by exotic long grasses {Pennisetum 
purpureum Schumach) and bamboos {B. 
tulda, B. balcooa, and M baccifera) for graz- 
ing of bison. An artificial lake and cashewnut 
(Anacardium occidentale L.) plantation on a 
small scale were other features of this site. 

Butterfly census 

Field surveys for butterfly fauna were con- 
ducted from May 2010 to October 2010 by a 
modified Pollard Walk Method (Pollard 
1977). One km long and 5 m wide permanent 
belt transects, one each in four sampling sites, 
representing a gradient of vegetation complex 
in four habitat types, were laid, and butterfly 
censuses were made between 08:00 and 11:00 
local time on four successive days. This was 
repeated at 30-day intervals, maintaining the 
same spatial scale in the four sampling sites. 
Thus, butterfly censuses were performed six 
times in as many months at each of the four 
habitats of TWLS. The cumulative area of 
monthly survey was 20 km^, which represents 
10.27% of the TWLS (194.71 km^). Point 
counts at 100-m-intervals were made along 
each transect to record butterfly species and 
their number. Walking pace was kept slow but 
at uniform speed, with a view to record the 
maximum richness and diversity. Stops were 



Table I. Transect information of four habitat types of Trishna Wildlife Sanctuary. 



Habitat type 


No. of transects 


Point counts 


Latitude 


Longitude 


Attitude (m) 


TWI.S I 


6 


60 


23''16'4L23" 


9r24'08.40" 


41 


TWLS II 


6 


60 


23"15'0L17" 


9r23'22.45" 


31 


I WLS III 


6 


60 


23"16'4L45" 


91"22'32.93" 


56 


TWLS IV 


6 


60 


23"16'33.29" 


91"23'01.29" 


52 
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Table 2. Number of individuals and species of butterflies recorded in Trishna Wildlife Sanctuary. 



Transects information 


Sampling sites 


TWLS 1 


TWLS II 


TWLS III 


TWLS IV 


No. of individuals 


298 


269 


264 


174 


No. of species 


39 


37 


33 


24 


Mean (± sem) no. of species/transect 


23.00 ±2.61 


22.17± 1.70 


19.33 ± 1.31 


15.33 ± 1.26 


Mean (± sem) no. of individual.s/transect 


49.67 ± 8.44 


44.83 ± 4.45 


44.00 ± 3.56 


29.00 i 4.23 



made along transect lines when an exact iden- 
tification of the species was not possible, and 
the specimens were then either photographed 
or captured by butterfly nets for closer exami- 
nation and identification. A few of the 
specimens that could not be identified in the 
field were brought to the laboratory for identi- 
fication, or were referred to the subject 
experts. During field study, collection data 
comprised of dates, location^ time, species, 
and number of individuals of each butterfly 
species were recorded separately for each 
sampling site. Vouchers of all field-collected 
specimens were deposited in the Ecology and 
Biodiversity Laboratories, Department of Zo- 
ology, at Tripura University. 

Butterfly identiflcation and geographic 
range 

Identification, nomenclature, and geographic 
distributions of collected butterflies followed 
Evans (1932), Wynter-Blyth (1957), Arora 
and Mondal (1981), Tsukada (1982), Otsuka 
(1988), Haribal (1992), Antram (2002), and 
Kehimkar (2008). Classifications of butterflies 
in this study were based on Ackery (1984). 
The geographic distribution ranges of differ- 
ent butterfly species were categorized on a 
numerical scale of 1-5 (smallest to largest) 
with some modifications to fit within the 
study area: (1) Eastern Himalayas (Sikkim to 
Assam); (2) Northeastern India and northern 
Indochina (up to Northern Myanmar); (3) In- 
do-Malayan region (India, including Andaman 
Island, Myanmar, and Malaysia); (4) Indo- 
Australian (Australasian tropics) including all 
India, Myanmar, and up to Srilanka; (5) 
Paleotropics (most of Asia north of the Hima- 



layan ranges apart from Africa north of the 
Sahara and all of Europe). 

Data analysis 

Raw data from the field were used to reveal 
species richness (Menhinik index), species 
diversity (Shannon- Weiner index), component 
of dominance (Simpson dominance index), 
and relative abundance of different species in 
a sampling site (Pielou's evenness index) 
(Magurran 1988). Comparisons in species 
composition between different forest habitats 
were estimated using single linkage cluster 
analysis based on Bray-Curtis similarity 
(McAleece 1998). The species recorded in this 
study were ranked based on relative abun- 
dance of individual species and also according 
to their known geographical range of distribu- 
tion. The raw data of species richness counts 
of six months from each sampling site were 
pooled to get rarefaction curves for compari- 
son of estimated species richness between the 
habitats. Sampling completeness was calculat- 
ed as a ratio of observed species richness to 
the richness estimate, and was expressed as a 
percentage (Sorensen et al. 2002). Biodiversi- 
ty Pro version 2 (Lambshead et al. 1997) was 
used to determine diversity indices, cluster 
analysis, species accumulation curve, rarefac- 
tion curves, and species richness estimates. 
Rank abundance diagram and species richness 
index were determined by Origin version 5 
(Microcal Software, Inc., www.microcal.com ) 
and PAST version 2 (Hammer et al. 2001), 
respectively. The Kruskal-Wallis test was 
used to compare species composition between 
different habitat types. 
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^able 3. List of butterfly species, their abundance, ranking, and geographic distribution ranking recorded in Trishna Wildlife Sanctuary. ^ 



Species 


F uinily 


Abundance 


Cico^^raphic ilistribiition ninkin^ 


1 WI.SI 


1 WISH 


rwi.siii 


rwLsiv 


Intal 


Arhopala pseudocentaurus (Doublcda> ) 


Lycaenidae 


2^ 


28 


18 


0 


68 


4 


Me Ian it is lec/a (L.) 


Nymphalidae 


16 


13 


21 


9 


59 


4 


Junonia almana (I..) 


Nymphalidae 


9 


9 


12 


18 


48 


4 


Orsotrioena medus V. * 


Nymphalidae 


15 


16 


10 


6 


47 


3 


Ypthima baldus 


Nymphalidae 


16 


10 


16 


4 


46 


3 


Catopsilia pomona (F.) 


Pieridae 


8 


9 


13 


14 


44 


5 


Euploea core (Cramer) 


Nymphalidae 


6 


8 


19 


10 


43 


3 


Scplis hy/as (1.) 


Nymphalidae 


5 


10 


15 


12 


42 


4 


Jitnonia allHes (L.) 


Nymphalidae 


18 


11 


0 


13 


42 


4 


Papilio polytes L. 


Papilionidae 


8 


10 


10 


11 


39 


3 


Ypthima /luehneri Kh by 


Nymphalidae 


9 


12 


12 


6 


39 


3 


Casta 11 us rosimon (F.) ^ 


Lycaenidae 


11 


9 


15 


0 


35 


4 


Narathwa selta (deNiceville) ^ 


Lycaenidae 


20 


13 


0 


0 


33 


3 


Catopsilia py rant he (L.) ' 


Pieridae 


7 


7 


9 


4 


27 


5 


Zizeeria karsandra (Moore) 


Lycaenidae 


5 


6 


8 


7 


26 


4 


Eiirema hccahe{L.) 


Pieridae 


4 


8 


8 


6 


26 


4 


Papilio de mole us F. 


Papilionidae 


12 


7 


0 


6 


25 


4 


lidoris Jarri (Moore) ^ 


Hesperiidae 


11 


9 


1 


0 


21 


4 


Pantoporia hordonia (stoll) ^ 


Nymphalidae 


4 


6 


5 


5 


20 




Lethe europa (F\) 


Nymphalidae 


8 


0 


11 


0 


19 


3 


Junonia hierta {¥.) 


Nymphalidae 


10 


0 


0 


8 


18 


3 


Zemeros flegvas ( Cramer) 


Lycaenidae 


7 


10 


0 


0 


17 


3 


Alraphancuru arislolochiae (F.) 


Papilionidae 


6 


6 


0 


5 


17 


3 


Tanaecia Icpidca (Butler) 


Nymphalidae 


4 


7 


4 


0 


15 


2 


Papilio nephelus Boisduval 


Papilionidae 


3 


0 


5 


7 


15 


4 


Zizina otis (F.) 


Lycaenidae 


8 


0 


0 


5 


13 


4 


Danaus genutia (Cramer) 


Nymphalidae 


3 


9 


0 


0 


12 


4 


A thy ma peri us (F.) 


Nymphalidae 


0 


4 


7 


0 


11 


3 


Mycalesis perse us (F.) 


Nymphalidae 


0 


-> 


5 


3 


10 


3 


Megisha malaya (Ilorsfield) 


Lycaenidae 


0 


4 


5 


0 


9 


3 


lamhrix salsala (Moore)* 


Hesperiidae 


9 


0 


0 


0 


9 


4 


SarangL'sa dasahara Moore 


Hesperiidae 


5 


0 


3 


0 


8 


4 


L daspes folus (C ramer ) 


Hespenidae 


0 


4 


0 


3 


7 


4 


Graph ium sarpedon (L.) 


Papilionidae 


0 


3 


4 


0 


7 


4 


Gandaca harina (Ilorstleld)* 


Pieridae 


7 


0 


0 


0 


7 


3 


Graphium agamemnon (L.)* 


Papilionidae 


0 


0 


7 


0 


7 


4 


Spalgis epius (W'cslw ood)* 


Lycaenidae 


0 


0 


0 


7 


7 


4 


Elymnias hypermnestra (L.) 


Nymphalidae 


3 


3 


0 


0 


6 


-) 


Troides helena (F.)*' 


Papilionidae 


5 


0 


0 


0 


5 


3 


Eurema hlamla (Boisduval) 


Pieridae 


2 


0 


0 


3 


5 


3 


Loxura atymnus (Sloll) * 


Lycaenidae 


0 


0 


5 


0 


5 


4 


Parnara gultatus (Bremer & Cirey) 


Hesperiidae 


0 


1 


2 


2 


5 


4 


Ilypolimnas holina (F.)* 


Nymphalidae 


0 


4 


0 


0 


4 


5 


Cepora nerissa (F.)* ' 


Pieridae 


0 


4 


0 


0 


4 


3 


Eraota timoleon (Stoll)* 


Lycaenidae 


0 


0 


4 


0 


4 


-) 


Arhopala Julia (Hewitson)* ^ 


Lycaenidae 


0 


0 


4 


0 


4 


-) 


Tagiades zapetus (Stoll) 


Hesperiidae 


-> 


-> 


0 


0 


4 


4 


Ancistroides nigrita (Fatrcillc)* 


Hesperiidae 


3 


0 


0 


0 


3 


2 


Lebadea martha (F.)* 


Nymphalidae 


3 


0 


0 


0 


3 


1 


/I riadnae ariadnae ( F . ) * 


Nymphalidae 


2 


0 


0 


0 


2 


3 


Discophora limora W'estwood* 


Nymphalidae 


0 


0 


2 


0 


2 


3 


Leplosia nina (F\)* 


Pieridae 


0 


0 


2 


0 


2 


4 


Discolampa ethion (Weslwood) 


Lycaenidae 


1 


1 


0 


0 


2 


4 


Elymnias malelas (Hewitson)*^ 


Nymphalidae 


0 


2 


0 


0 




2 


Papilio memnon L* 


Papilionidae 


1 


0 


0 


0 




3 


Euhalia aconthea (Cramer)* " 


Nymphalidae 


0 


1 


0 


0 




3 


f.ethe me car a (Moore)* 


Nymphalidae 


0 


0 


1 


0 




"> 


Catapaecilma elegans (Druce)* 


Lycaenidae 


0 


0 


1 


0 




4 


Suasfus gremius (V.)* ' 


Hesperiidae 


0 


1 


0 


0 




4 



Signs with species names denote the following: ^unique species, - singleton species, + present in all sampling sites, x threatened species. 

V 



Results 



Butterfly fauna 



A total of 1005 butterflies representing 59 
species belonging to 48 genera and 5 families 
were recorded in this study (Table 2, 3). A 
maximum of 39 species and 298 individuals 
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Table 4. Diversity parameters and species richness estimates of butterfly communities in the four habitat types. 


Habitat t>-pes 


Menhinik index (3//!) 


Shannon-W einer diversitj index (//,) 


Simpson's dominance index (Z>.) 


Pielou's evenness index (J') 


Chao2 


Observed species richness 


Sampling completeness (%) 


TV\ LS 1 


2.26 


1 .50 


0.03 


0.94 


39 


39 


100 


TWLS II 


2.26 


1.47 


0.04 


0.94 


39 


37 


95 


TWLS III 


2.03 


1.41 


0.04 


0.93 


34 




97 


TWLS IV 


1 X2 


1 32 


005 


096 


24 


24 


100 


Oveiall 


1.86 


3.62 


0.03 


0,89 


61 


59 


97 



of butterflies were recorded in TWLS-I, which 
was dominated by mature mixed moist decid- 
uous forest, and minimum of 24 species and 
174 individuals were recorded in TWLS-IV, 
which was dominated by long grasses and 
bamboos (Table 2). Sampling sites TWLS II 
and TWLS III were represented by 37 and 33 
species, respectively, but without any signifi- 
cant difference (Kruskal-Wallis Test: H = 
0.18, p = 0.681) between the two sites. 
Among the four sites, TWLS I and TWLS II 
showed significant differences with TWLS IV 
in terms of species composition (TWLS I: H = 
6.19, p = 0.018; TWLS II: H = 4.52, p = 
0.045). Out of 59 species, 21 species (35.59%) 
were recorded from one of the four sampling 
sites only and, hence, these were considered 
'unique' species for the purpose of this study. 
These included five 'singleton' species (Table 
3). Another 14 species (23.73%) were found 
to be common to all the sampling sites. Nine 
species (15.25%) were recorded as 'threat- 
ened' or 'endangered' as per Part IV of 
Schedule I and Schedule II of the Indian 
Wildlife Protection Act 1972 (Anonymous 
2006), and another one species, Troides Hele- 
na (L.), was listed under CITES Appendix II 
(Collins and Morris 1985). 

Butterfly diversity 

Family- wise distribution of butterflies showed 
that members of Nymphalidae dominated the 
collection (23 species and 492 individuals) 
followed by Lycaenidae, Papilionidae, Hespe- 
riidae, and Pieridae, in that order (Table 3). 

The diversity parameters of butterflies showed 
variations in the four sampling sites. In gen- 
eral, the four sampling sites showed high 
species richness and diversity of butterflies. 



very low dominance of species, and high 
evenness of distribution (Table 4). Specifical- 
ly, the mature mixed moist deciduous forest at 
TWLS I showed maximum diversity {Hs= 
1.50) and least dominance (A = 0.03) of but- 
terflies, whereas the habitat with long grasses, 
bamboos, and plantation crops (TWLS IV) 
showed minimum diversity {Hs= 1.32) and 
highest dominance (A = 0.05) of butterflies. 
Evenness of distribution in all the study sites 
was found to be high {J= 0.93 to 0.96). Spe- 
cies diversity in regenerated mixed moist 
deciduous forest with little understory vegeta- 
tion (TWLS II) and in secondary mixed moist 
forest with patches of bamboos (TWLS III) 
were lower than TWLS I but higher than 
TWLS IV. 

Species ranking 

A ranking of 59 species according to their ge- 
ographical distribution showed that 26 species 
were distributed in the Indo-Malayan region 
(rank 3), another 22 species in the Australa- 
sian tropics (rank 4), 3 species in Paleotropics 
(rank 5), and the remaining 8 species showed 
restricted geographical distribution, which in- 
cluded 7 species distributed in northeastern 
India and northern Indo-China up to northern 
Myanmar (rank 2), and the other species, La- 
badea martha (F.), was distributed in the 
eastern Himalayas (rank 1) only (Table 3). 

The distribution of ranks according to the 
abundance of different butterfly species 
showed that 59 species were distributed into 
33 ranks. Amblypodia centaurus (Doubleday) 
and Melanitis leda (L.) showed relative abun- 
dances of 6.77% and 5.87% respectively. 
Together, these species accounted for 127 
(12.64%) of the individuals encountered in 
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this study. The next 21 species showed abun- 
dances in the range of 1.69-4.78%, and 
together accounted for 67.06% of the total 
abundance. Another 31 butterfly species to- 
gether accounted for 19.80% of total 
abundance. The remaining 5 species were 
found to be singletons (Table 3), and these 
contributed 0.05% of the total abundance 
(Figure 2). 

Richness estimates 

Species richness estimate, using Chao 2, was 
found to give the best estimate for the samples 
of this study. Estimation of species richness in 
the four habitats of TWLS by Chao 2 showed 
expected richnesses that were very close to the 
observed values (Table 4). The overall esti- 
mate of species richness from the four habitats 
was higher by 2 species than the observed 
value. This was also evident from the high 
values of the sampling completeness of this 
study, which varied in a narrow range be- 
tween 95% and 100%) between four habitats, 
and which was 97% overall. Rarefaction 
curves from the four habitats showed quick 
rises at first and than either leveled off 
(TWLS IV) or approached asymptote gently 
(TWLS I, II, and III) (Figure 3). 

Butterfly species composition 

The cluster analysis based on Bray-Curtis sin- 
gle linkage similarity index demonstrated the 
differences and similarities between the but- 
terflies species composition recorded in four 
habitat types. The open habitat of long grasses 
and bamboos (TWLS IV) stood out clearly 
from the other 3 habitats and showed linkage 
at 50.86% similarity, which represented the 
lowest similarity. The degraded forest habitat 
(TWLS III) was linked at 58.16% similarity to 
the cluster of habitats of primary forest 
(TWLS I) and regenerated forest (TWLS II), 
which showed the highest similarity in species 
composition (70.55%) (Figure 4). 
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Discussion 

This is the first study on the distribution and 
abundance of butterflies in the quasi- 
evergreen and moist deciduous forests of 
TWLS. The 59 recorded species in this sanc- 
tuary compares favorably with the 71 species 
found in Aralam wildlife sanctuary in Kerala, 
south India, which has similar habitats 
(Sreekumar and Balakrishnan 2001). All the 
species recorded in Aralam wildlife sanctuary 
(55 sq. km) were 'generalist' and none of 
them were considered to be threatened or of 
rare occurrence. In contrast, eight species, i.e., 
Lethe europa (F.), Cepora nerissa (F.), Casta- 
lius rosimon (F.), Narathura selta (de 
Niceville), Pantoporia hordonia (stoll), Me- 
gisba malaya (Horsfield), Arhopala fulla 
(Hewitson), and Baoris farri (Moore) record- 
ed in TWLS were enlisted in the Indian Forest 
Act 1972 under Schedule I & II, and another 
species, Troides Helena L., was listed in Ap- 
pendix -II of CITES, and these are protected 
by law. 

Overall species abundance and richness re- 
vealed that Nymphalidae was the most 
speciose and individualized family, Pieridae 
was the least speciose family, and Hesperidae 
was the least individualized family in the 
study area. In the context of tropical environ- 
ment, similar patterns of species abundance 
and richness were reported from the western 
ghats and western coast of India (Kunte 1997; 
Eswaran and Pramod 2005; Padhye et al. 
2006; Krishnakumar et al. 2007; Raut and 
Pendharkar 2010), from the Silent Valley Na- 
tional Park, Kerala, south India (Mathew and 
Rahamathulla 1993), and in the Parambikulam 
wildlife sanctuary, south India 
(Sudheendrakumar et al. 2000). The domi- 
nance of Nymphalidae species may be 
attributed to their being polyphagous, which 
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helps these butterflies to live in a variety of 
habitats, and also because many species of 
this family are active fliers, which helps them 
to forage larger areas. However, in the pre- 
sent study, Pieridae was poorly represented in 
comparison to studies made in other sanctuar- 
ies of India, which reported members of 
Hesperiidae to be least represented. One of the 
possible reasons for this difference could be 
due to the difficulties in observing Hespe- 
riidae butterflies because of their dull color 
and ability to fly rapidly following any dis- 
turbance. 

Fifty-nine species recorded in this study 
showed unequal distribution of abundance in 
the four habitats of TWLS. Several studies 
have reported the influence of habitat disturb- 
ance on abundance of butterflies (Hill et al. 
1995). In the present study, 13 species having 
a geographical distribution ranking of 3-5 
were present in all four habitats. Another 21 
species showed low abundance and habitat 
speciflcity (Table 3). TWLS IV contained one 
of these, Spalgis epius (Westwood), whereas 
TWLS I, TWLS II, and TWLS III contained 
7, 5, and 8 unique species, respectively. How- 
ever, 6 of these were truly restricted in their 
geographic distribution, and the other 15 spe- 
cies showed wide geographical distribution. 
Geographically restricted species are consid- 
ered habitat speciflc, show low ecological 
tolerance, prefer undisturbed to least disturbed 
habitats, and have high conservation value in 
comparison to widely distributed species, 
which show high ecological tolerance, occur 
in a gradient of vegetation complex, share 
food resources over a wide geographical 
range, and are of low conservation value 
(Spitzer et al. 1997). In the present study, spe- 
cies poor habitat TWLS IV contained one 
unique species, S. epius, which is aphi- 
dophagous. The results clearly showed that 
both species rich habitats and species poor 
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habitats can support unique butterfly species 
depending on their resource utilization. 

Habitat selection in butterflies is directly re- 
lated to the availability of preferred food 
plants for larvae and adults (Grossmueller and 
Lederhouse 1987; Thomas 1995). In the pre- 
sent study, the maximum number of species 
and individuals were observed in mature 
mixed moist deciduous forest (TWLS I), fol- 
lowed by regenerated secondary mixed moist 
deciduous forest (TWLS II). 

The degraded forest of TWLS III, with patch- 
es of woods, grasses, and herbs, supported 33 
species of butterflies. The presence of plenty 
of woods and dense understory delimited by 
profuse growth of diverse herbs and shrubs in 
these habitats provide a rich nectar source to 
adult butterflies and food to larvae. The poor 
species richness and low diversity of butter- 
flies recorded in the grassland habitats of 
TWLS IV corroborate with the flndings of 
Ramesh et al. (2010), who recorded a similar 
distribution pattern of butterflies in the west- 
em coasts of India. The cultivation of exotic 
grass species, which creates problems for the 
host-speciflc butterfly species that rely on lo- 
cally available plant species for their survival, 
may be responsible for the poor species rich- 
ness and diversity of TWS IV. Similarly, 
invasive species impacted many native herbi- 
vores, particularly those species for which the 
plant serves as a potential food plant (Nagy et 
al. 1998). The high species richness, abun- 
dance, diversity, and uniqueness of TWLS I 
make it a key habitat for TWLS. 

The results showed that more than 50% of the 
butterfly species recorded in all the habitat 
types was the same, despite differences in 
habitat features. The high similarity values of 
butterfly fauna between different habitats of 
TWLS indicate low beta diversity over short 
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distances between the four transects in this 
forest, and bear resemblance to a study con- 
ducted in tropical forests, where herbivorous 
insects showed low beta diversity even over 
large areas (Novotny et al. 2007). It is gener- 
ally accepted that rich floral diversity in 
tropical forests promotes herbivores, many of 
which are generalists (Price 1997). 

Many studies have reported the relationship 
between habitat heterogeneity and species di- 
versity (Bazzaz 1975; Brooks 1997; Atauri 
and Lucio 2001; Tews et al. 2004). In most 
habitats, plant communities determine the 
physical structure of the environment, and 
therefore have considerable influence on the 
distributions and interactions of animal spe- 
cies (Lawton 1983; McCoy and Bell 1991). 
The highest diversity, minimum dominance, 
and occurrence of seven unique species of 
butterflies recorded in TWLS I may be at- 
tributable to its vegetation complexity and 
multilayered canopy, which facilitate different 
sets of microclimates, making the habitat dis- 
tinct for different butterfly species. 
Rosenzweig (1981) stated that diversity is en- 
hanced by the presence of specialists that 
exhibit distinct habitat preferences. TWLS II 
also attracted 37 species of butterflies due to 
the edge effects of the neighboring habitats of 
TWLS I and TWLS III, and an additional fea- 
ture of a large number of pits caused by 
uprooting of plants by wild boars. Pits provide 
sites for deposition of mineral rich feces and 
urine perfect for mudpuddler butterflies, such 
as Papilio polytes L., Papilio demoleus L., 
Castalius rosimon Fruhstofer, and Eurema 
hecabe (Moore). This flnding is in agreement 
with the flndings of Ramos (1996). 

The species estimator Chao 2 order was found 
to be a reliable overall predictor with respect 
to accuracy of species richness in the four 
habitat types of this study. Peterson et al. 
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(2003) and Brose et al. (2003) reported similar 
results using Jackknife 2 order for small sam- 
ple sizes. 

The flndings indicated that mature secondary 
and regenerated forests supported high butter- 
fly diversity and species richness, while exotic 
grassland had a negative impact on butterfly 
community. The occurrence of six threatened 
species and 21 unique species in TWLS 
strongly suggests that there is a need for the 
implementation of sustainable conservation 
strategies for the protection of the rare taxa of 
the sanctuary. 
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Figure I. Location map of Trishna Wildlife Sanctuary, High 
quality figures are available online. 
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Figure 2. Rank abundance diagram of butterfly species record- 
ed at Trishna Wildlife Sanctuary. High quality figures are available 
online. 




Figure 3. Sample based rarefaction curves of estimated species 
richness at four sampling sites of Trishna Wildlife Sanctuary. 
High quality figures are available online. 
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Figure 4, Single linkage cluster analysis between habitat types 
based on Bray-Curtis similarity; abbreviations denote: TWLSI = 
Trishna Wildlife Sanctuary I, TWLSI! = Trishna Wildlife Sanctu- 
ary II, TWLSIII = Trishna Wildlife Sanctuary III, TWLSIV = 
Trishna Wildlife Sanctuary IV. Dissimilar letters following habitat 
types indicate significant differences by Kruskal-Wallis test {p < 
0.05). High quality figures are available online. 
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